FORCE INDUCED DEPINNING OF DIRECTED POLYMERS 

GIAMBATTISTA GIACOMIN AND FABIO LUCIO TONINELLI 

Abstract. We present an approach to studying directed polymers in interaction with 
a defect line and subject to a force, which pulls them away from the line. We consider in 
particular the case of inhomogeneous interactions. We first give a formula relating the 
free energy of these models to the free energy of the corresponding ones in which the 
force is switched off. We then show how to detect the presence of a re-entrant transi- 
tion without fully solving the model. We discuss some models in detail and show that 
inhomogeneous interaction, e.g. disordered interactions, may induce the re-entrance phe- 
nomenon. 
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1. Introduction 

1.1. Overview. Pulling polymers out of potential wells by applying a force / in the 
direction orthogonal to the pinning region has been considered at several instances (see 
e.g. |111 1131 1141 115j and references therein). The techniques employed essentially lead 
to exact (and mathematically rigorous) results for models with homogeneous interactions 
(see in particular |14j and |15j). while in the inhomogeneous case, notably in the disordered 
case, the results are often based on arguments out of mathematical control (e.g., replica 
and renormalization group computations) or on annealed approximations that, in general, 
are not close to the behavior of the quenched system. One of the interesting phenomena 
that has been pointed out, at least in some model systems, is the presence of a re-entrant 
transition (see in particular |14j . but this issue is taken up in most of the references we 
have given). By this we mean that, for a fixed force /, the polymer is pulled out of the 
defect line at low temperature and at high temperature, but for intermediate temperatures 
it is localized at the defect line, very much like if the force was not present. References 
above are limited to theoretical work, but all of them are motivated by real experiments: 
the literature on real experiments is very rich and here we single out that deals with 
DNA unzipping and contains further references. 

Our purpose is to point out that one can directly relate the free energy in presence 
of a force to the free energy of the polymer with free endpoint for a very large class of 
directed models and that the presence of a re-entrant transition is easily related to suitable 
asymptotic behavior of the expressions appearing in this formula. As a consequence, we 
will give a simple necessary ad sufficient condition for the re-entrance to take place for 
very general models (including disordered models). 

In order to be concrete we choose to deal with a precise and rather limited class of 
models, but the reader will realize that the method we present is very general. Extensions 
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are considered explicitly in Section |5J Moreover, the phenomenology of the restricted class 
of models we consider is already very rich. 

We point out that we focus exclusively on cases in which the polymer interacts only 
at the defect line and other relevant cases, like the case in which bulk disorder is present 
(e.g. jjj), are not considered here. On the other hand, it is straightforward to generalize 
the content of this note to directed models of copolymers near a selective interface [HI E3 • 

1.2. A model. Let us consider a (p, g)-walk, that is a random walk S := {S n } n= o,i,... with 
independent identically distributed (IID) increments {S n+ \ — 5' n } ri= o,i,..., 5*0 := and Si 
taking values in {-1,0, +1}, with P(Si = 1) = P(5i = -1) = p/2 and P(5i = 0) = q. 
We assume p+q = 1 and p > 0. Particularly relevant for what follows is the distribution of 
r := inf{n : S n = 0}: therefore we set K(n) = P(r = n), along with K{n) := Ylj>n^U)- 
We point out that the Laplace transform of K(-), that is ^ n exp(— bn)K(n), b > 0, can 
be written explicitly and K(n) itself can be expressed in terms of P(Sk = 0), k = 1, . . . , n, 
by using a standard renewal theory formula 

n 

P(S n = 0) = 6 n , + K(j)P(S n -j = 0). 

In particular one obtains that n 3 / 2 K(n) converges to c p := vV( 27 8, App. A], if p < 1 
(an analogous result holds for the p = 1 case, see below). Note also that K(0) = 1, since 
S is recurrent. 

Given a sequence u> := {uj r i} n of real numbers (the charges) we define the new measure 
Pjy w on the space of random walk paths by introducing a Boltzmann weight {(3 > 0, 
/ > 0): 

d p/3J / N \ 

—^(S) = -— exp /?;>> n l 5n= o + /3fS N . (1.1) 

We will consider two (wide) classes of charges: 

(1) w is a periodic (deterministic) sequence, i.e. a sequence of numbers such that 
^n+T = w n for some positive integer T and every n. We denote by t(w) the minimal 
T with such a property, that is the period. If t(lo) = 1 then oj is homogeneous. 

(2) oj is the realization of a sequence of random variables. For simplicity we consider 
only the case of IID variables, therefore the law of oj is determined by the law of 
oji, but our analysis would go through in the much wider domain of stationary 
sequences of variables. In the end, the properties of the system are determined by 
the law P of oj. We refer to this case as disordered or quenched. 

The existence of the limit 

km ilog4° = : HP), (1.2) 

N—>oo IV 

is very well known. This follows from elementary super-additivity arguments in the pe- 
riodic (and, of course, in the homogeneous) set-up (see e.g. [HI Ch. 1]): in general, the 
limit depends on w. In the disordered case instead one has to be a bit more careful: a 
precise statement is for example that if M\oji\ < oo then the limit in (|1.2|) exists P-almost 
surely and also in the L 1 (P) sense. Moreover the right-hand side is in principle a random 
variable, but it turns out to be degenerate, that is almost surely independent of oj. This 



DEPINNING AND RE-ENTRANT TRANSITIONS 



3 



property, usually referred to as self-averaging, is well understood in this context (see jSJ 
Ch. 4] for proofs and overview of the literature). We remark that 

> E [S n + 0, forn = 1, 2, . . . , N] = ~K(N) > cN~ 1 ' 2 , (1.3) 

for some c > 0, so that f(/?) > for every (3. Much literature has been spent on this 
model and about the fact that f(/3) > corresponds to localized regime, that is to the 
case in which the typical trajectories are tightly bound to the defect line (cf., e.g., |10j 
and |H[ Ch. 7]). On the other hand f(/3) = corresponds to a delocalized regime. One 
should however distinguish between the critical point (3 = j3 c = sup{/3 : f(/3) = 0} and the 
truly delocalized region where (3 < (3 C . Several results are available also on the delocalized 
regime (see [H] and references therein, in particular |Hj for the disordered case) and they 
say, very roughly, that very few visits are paid by the polymer to the defect line. 

A relevant difference between the two classes of charges we consider is that in the first 
case f(/3) is explicitly known. To be more precise in the periodic case f(/?) can be expressed 
in terms of the leading eigenvalue of a suitable t(uj) x t(uj) Perron-Frobenius matrix [31 
(see also Appendix I A. 2|) : computing such an eigenvalue becomes harder and harder for 
larger periods, but the problem trivializes in the homogeneous case. In the disordered case 
instead only estimates on f(/3) are known. 

We point out also that for every the limit 

lim 1 logE [exp (yS N ) ; S n + forn = 1, 2, . . . , N] := G (y), (1.4) 

TV^oo A 

exists. Actually, the explicit value of g(-) is easily computed by applying the reflection 
principle [SJ App. A] and one obtains G(y) = log(pcosh(y) + q). 
For the model with the force we have: 



Proposition 1.1. For every (3 and f the limit, that we denote by f(/3, /), of the sequence 
|(1/A) log Z^-'^l exists (W(du)- almost surely (a.s.), in the disordered case). Moreover 
we have the formula 

F(P,f) = max(F(/3),G(/3/)). (1.5) 

A line of non- analytic points of f(-, •) is therefore evident: for every (3 we set f c {(3) '■= 
f3~ 1 G~ 1 (f(/3)) G [0,oo) (we are looking at g(-) as a function, and a bijection, from [0,oo) 
to [0, oo)). Therefore for every (3 

mf) = \cW) if/>«/3). (L6) 

It is quite easy to get convinced that this non-analyticity corresponds to a localization- 
delocalization transition: the most interesting case is when f(/3) > 0, otherwise the system 
is already delocalized at / = 0. By convexity, <9/-f(/3, /) > if / > f c {(3) (and if dfF{(3, f) 
exists). This directly implies that lini/v^oo Ei^^Sn/N > (and the same statement holds 
even if df¥{(3, f) does not exist, but one has to replace the lim with liminf). For / < f c {0) 
instead lim7v-*oo E^^Sat/A = 0. These are distinctive marks respectively of localization 
and derealization (for sharper results we refer to |16j). 

It is also worthwhile to observe that formula (|1,5[) yields that df¥{(3, f) is discontinuous 
at / = f c (f3) for (3 > 0. So the transition is of first order (as argued in much of the previous 
literature). The underlying mechanism has also been expoited in |2j. 



1 



GIAMBATTISTA GIACOMIN AND FABIO LUCIO TONINELLI 



We now turn to the issue of the existence of a re-entrant transition. As mentioned above, 
this refers to the fact that for some fixed / the system undergoes (at least) two phase 
transitions as the temperature T = 1/(3 is increased from zero to infinity. More precisely, 
two cases are observed: one can either observe a pattern of the type localized-delocalized- 
localized by increasing T, or the opposite one: delocalized-localized-delocalized. 

Of course, re-entrance is equivalent to the non-monotonicity of f c as a function of (3 and 
a sufficient condition for it is that the difference / c (oo) — / c (0) has the opposite sign as 
9/3 / c (0). As we show in the following, these quantities are related to the asymptotic be- 
havior of f(/3) for (3 — > and (3 — > oo, which in many cases (including quenched disordered 
situations) can be easily computed without fully solving the model. 

We wish to emphasize that, while in principle the occurrence of 

^/c(0)[/ c (oo)-/ c (0)]<0, (1.7) 

is just a sufficient condition for re-entrance, it turns out numerically in the cases we have 
checked (cf. Fig. ^and|2j) that when ljl.7JI is not verified the critical force is monotone in 
(3 and re-entrance is absent. 



2. The homogeneous case 
This section has the two aims: 

(1) generalizing the results of the previous literature. For example in ^S] only the 
cases q = 1/3 and q = are considered and no re-entrance is observed, while for 
different values of q re-entrance does appear; 

(2) providing the general scheme that we follow also in the next sections, even if the 
models of this section are exactly solvable (but we will not solve them). 

First of all note that, if uj n = — c < the polymer is delocalized already in absence of 
the force, and f c {(3) = for every [3 > 0, an uninteresting situation. Therefore, in the 
homogeneous case, we will assume that uj n = c > and, without loss of generality, c = 1. 

The behaviors of the critical force for (3 \ and for (3 — * oo are easily obtained: 
Proposition 2.1. For the homogeneous model, i.e. Ql.ljl and uj n = 1, with p < 1 one has 

Hm/„<ffl = \, (2.1) 

and 

fciP) =°°l + ^ (log? - logb/2)) + o (J) . (2.2) 

In particular, one has re- entrance if and only if 2/3 < p < 1. As for p = 1, one has 
f c {(3) -> 1 for f3 -> and 

A(«^i + i.o g 2 + (i). (2.3) 

Re-entrant behavior is not observed in this case. 

See also Figure ^a), were we plot f c {(3) as a function of 1/(3, for different choices of q. 
Note that the transition p — * 1 to the simple random walk case is singular. 
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Figure 1. Critical force f c as function of 1//3, for the homogeneous model, (a) Without 
hard wall. The curves correspond q — 1/2, 1/3, 1/5, 1/50 and 0, from top to bottom. 
Re-entrance is present if and only if < q < 1/3. The critical force / c (/3) is com- 
puted by using 11.51 . where f(/3) is the solution of K(n) exp(— F(/3)n) = exp(~/3) 
|5] Proposition 1.1]. (b) With hard wall (the values of q are the same as in (a)). Re- 
entrance is present if and only if q > 1/3. In this case, f(/3) is given by the solution of 
E n ^ + (n)exp(-F(/3)n) = exp(-/3), where K+(n) = K{n)/2 if n > 2 and K+(l) = 
K{1). 

Proof of Proposition \2.1\ Consider first the case p < 1. Since K(n) ~ c p n -3 / 2 , c p = 
\Jpj (2ir), for \ one has |HJ Theorem 2.1] 

F (0) r (2 4) 

HP) (2c p r(l/2))2- w 

This, together with a Taylor expansion of G(y) around y = 0, immediately gives ()2.1j) . 

The large behavior can also be easily captured without explicit computations: for 
p < 1, f(0) = + logK(l) + o(l) as — > oo, and K(l) = q = 1 — p. This corresponds 
simply to the fact that for large the dominant trajectories are the ones such that |{n : 
S n ^0}\ <. N. Since G(y) = y + log(p/2) + o(l) as y -> oo, one finds (j^ . 

For p = 1, (|2.4|) is modified in this case into f(0) ~ 2 /2 (this follows from B, Theorem 
2.1] plus the fact that, for the simple random walk, K(2n) ~ n - 3 / 2 v /1/(4tt) Ch. Ill] 
and K(2n + l) = 0). As a consequence, f c (P) — > 1 for /3 ^ 0. As for the — > oo behavior 
notice that, f(/5) = 0/2 + (log#(2))/2 + o(l) and iT(2) = 1/2, so that (E3J) follows. 
The statements about re-entrance easily follow from ()2. 1(1 - 1)2.31) . □ 

2.1. The model with hard wall repulsion. We conclude this section by showing that 
the phase diagram and the re-entrance phenomenon are strongly model dependent. To 
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this purpose, we modify the model by inserting a hard wall condition (like in ^5])> which 
corresponds to inserting in the right-hand side of the indicator function of the event 
{S : S n > for n = 1, . . . , N}. Of course, the partition function and the free energy will 
be in general modified. We consider for conciseness only the case p < 1. 

Proposition 2.2. For the homogeneous model with hard-wall repulsion and p < 1, there 
exists flo > such that f c (f3) = for (3 < (3q. Moreover, the large j3 behavior of f c ((3) is 
still given by ()2.2|) . Re-entrance takes place for p < 2/3. 

Proof of Proposition^^ In presence of the hard wall, for (3 sufficiently small but finite one 
has f(/3) = (cf. for instance [El Section 1.2]), so that / c (/3) = 0. This is rather intuitive: 
even in absence of the pulling force, the entropic repulsion effect provided by the wall is 
enough to delocalize the polymer. As for the [3 — ► oo limit, since the dominant trajectory 
S n = is allowed by the hard wall constraint, one finds again f(/3) = j3 + log K(l) + o(l), 
and Equation (|2.2|) . Occurrence of re-entrance for p < 2/3 immediately follows from the 
small- and large-/? behavior of the critical force. □ 
Observe that for the model with hard wall repulsion the situation is somewhat reversed 
with respect to the previous case: one has re-entrance for p < 2/3 and no re-entrance for 
2/3 < p < 1 (see also Fig. Ufb)). 

3. The periodic case 

We start by recalling that the solution of the periodic case can be reduced to a finite 
dimensional, in fact T(w)-dimensional, problem (0 [IJ [H]): by solving numerically this 
finite dimensional problem we have drawn the curves in Fig. [21 However, from the Tri- 
dimensional problem one can extract many analytic features too. Here we concentrate on 
small and large (3 behavior, since they suffice to highlight part of the variety of observed 
phenomena: we collect and discuss here the results, see App. lA~2l for the proofs. 

Just a bit of notation for the following result (for sake of conciseness, we restrict our 
attention to the case without hard-wall repulsion). Set £q := min{n : oj n = +1} (we are 
implicitly assuming that uj n is not equal to —1 for all n) and := inf{n > l{ : uj n = +1}. 
Set also i(lu) := min{i : t- L > t(w)}. 

Proposition 3.1. In the periodic case we have the small temperature behavior 



^V(2a, n -l) + i 4tE ] 



(3.1) 



feW ^ ^>](^n-l) + i[ ^>>g#(4-4_l)-]0g(p/2) | +0(1/(3) 

and at high temperatures (/3\,0) 

fc(P) 



^ ^ + o((3) if E?=1 ^ n > 0, 

kP + o{(3) i/ElL1^ = 0, (3.2) 



2p 

.0 i/En=i w n<0 and/3<p , 

where (3q := sup{/? : f((3) = 0} ; and f3 £ (0, oo) if En=i w n < 0. 

Of course such a result is sufficient in order to check our sufficient condition (|1.7j) for 
the occurrence of a re-entrant transition, for a given periodic charge sequence lo. On the 
other hand, a full characterization of the critical curve in general requires some numerical 
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computations (but we stress once again that they are just finite dimensional computa- 
tions). Note however, by comparing with the homogeneous case, that the inhomogeneous 
character of the charges induces a re-entrant transition for example if q = 1/3 (but of 
course it can induce it also for q < 1/3). 




Figure 2. Critical force f c as function of 1//3, in the case of periodic charge distribution 
and q = 1/3. In the three cases t(u) = 50 and the charges have been chosen at random, 
i.e. each sequence {^n}n=i,...,50, w„ 6 { — 1,1}, is equiprobable (the explicit sequences 
{(u> n + l)/2} n= i i ... j 5o are given next to the corresponding curves): for the top curve 
\{n : uj n = +1}| = 30, while for the bottom one \{n : uj n — +1}| = 22. The intermediate 
curve corresponds instead to the case in which there are as many positive charges as 
negative, and in fact the curve approaches for large temperatures. As explained in the 
text, these graphs give also strong hints on the behavior of disordered systems. 



An important feature resulting from Prop . |3~T1 (but also from Fig.[2J) is the strong depen- 
dence on oj. In particular the sign of the mean over a period leads to drastically different 
behaviors for large temperatures. It should be however pointed out that if u%, . . . , u r are 
sampled in a IID fashion, like in the quenched disordered case, then in the limit of T — > oo 
the free energy of the periodic model converges P(do;)-a.s. toward the free energy of 
the disordered model [HI Th. 4.5]. One then sees directly that such a result implies the 
convergence of the critical force of the model with period T to the critical force of the 
corresponding disordered model. In Fig. [21 we have plotted a case in which X^n=i w n = 
and two cases deviating above and below the mean by more than one standard deviation. 



4. The disordered case 



Let us now have a look at the disordered model. For the sake of conciseness, we consider 
only the case p < 1 and P(wi = +1) = P(u>i = -1) = 1/2. 
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It is known that in this case the model with / = and without hard wall repulsion 
is localized for every j3 > (see [HI Ch. 5]). The behavior of the free energy is not fully 
under control, but one can prove the following: 

Proposition 4.1. For the disordered model, f c (P) —> for (3 — > 0. As for large (3, one 
has 

m 1 ■ 1 ( 1 



\ + \ ^E2- J logK(j)-log(p/2)j +o(£) (4.1) 



so that re-entrance is observed as soon as 



-. oo 

-^V'logtftf) > log(p/2), (4.2) 



2 

i=i 

that is for q > 0.1994.... 

Proof of Proposition 14.11 For (5 — > 0, the annealed bound 

F^^^llogE^, (4.3) 

is already sufficient to prove that / c (/3) — > 0. In fact, the right-hand side is just the free 
energy of the homogeneous model where (5 is replaced by log(cosh /3), and therefore (recall 
Q2.4|) and discussion thereafter) it behaves like /3 4 / (8p) for /3 small. This immediately 
implies that f c (/3) -> for /3 0. 

Remark 4.2. In 1 it is actually proven that f(/3) p ^ /? 4 /(8p), so i/tai f c {(3) ~ 
/5/(2v^)- 

As for large /?, we need the following: 
Lemma 4.3. For every e > £/iere exists j3 £ such that for (3 > j3 e we have 



f(/3)-| --^2-n og ^(i) 



< £■ (4.4) 



Note that (|4,1|) is a direct consequence of Lemma 14.31 and the proof of the latter can 
be found in Appendix IA.11 Here, we just notice that ()4.4|) shows how the situation is 
different from the homogeneous case, where one has rather f(/3) = (3 + \ogK{\) + o(l). 
In fact also in the disordered situation the dominant configurations are those for which 
S n = for every n such that oj n = +1 and S n ^ otherwise. However, now the average 
density of +1 charges is 1/2 (whence the leading term /3/2) and the distance between two 
successive positive charges is random and can take any value j = 1,2,... with geometric 
probability 2~ 3 . □ 

As explained in Section |31 we expect the graph of f c (-) to be close to the intermediate 
curve in Fig. |2 



DEPINNING AND RE-ENTRANT TRANSITIONS 



9 



4.1. The model with repulsion. Like in the homogeneous case, if we add the hard-wall 
repulsion then for f3 sufficiently small f(/3) = f c (P) = (see for instance (HI Proposition 
5.1]). Moreover, a look at the proof of (|4.4|) shows that Lemma 14.31 still holds in this case, 
provided that K(n) is replaced by K + (n) defined as K + (l) := K(l) and K + (n) := K(n)/2 
for n > 2. This is just related to the fact that, if n > 2, of all the possible trajectories 
{So, . . . , S n } of the (p, q) walk satisfying So = S n = 0, S{ ^ for 1 < i < n, only half 
survive the introduction of the hard wall constraint. Therefore, (|4.1j) still holds with K{.) 
replaced by K + {.) and a sufficient condition for having a re-entrant transition is 

^ oo 

-^2-nog^+(i) > Iog(p/2). (4.5) 

i=i 

Numerically, this corresponds to q > 0.2838 .... 

5. Generalizations 

The procedure presented in this paper to study the phase diagram and the occurrence 
of a re-entrant phase transition applies well beyond the case of the (p, q)-walks. We have 
seen in fact that all what one needs to know is the asymptotic behavior of G(y) for small 
and large y, and the asymptotic behavior of f(/3) for small and large (3. 

Just to show an example of generalization, consider the case of a directed polymer 
in (2 + 1) dimensions with homogeneous {uj n = 1) pinning attraction to a defect line. 
In particular, we will assume that S n = (Sn\Sn^) € 1? and {S' n } n= o,i,... is the two- 
dimensional simple random walk where S n +\ is chosen uniformly among the 4 neighbors 
of S n . In this section, it will be understood that K(n) is the probability that the first return 
to zero of the two-dimensional random walk occurs at time n (in the previous sections the 
same symbol was used for the analogous quantity referring to the (p, q) -random walk). 

Again, the Boltzmann weight is defined as in Ql.l|) . with the only difference that the 

term fSjv is replaced by, say, fsff. This corresponds to assuming that the force is 
pulling in the direction "1" (orthogonal to the defect line). Then, Proposition 11.11 still 

holds, provided that in the definition of G(y) one puts yS$ instead of ySw- 
The analogue of Proposition 12.11 is the following: 

Proposition 5.1. For the homogeneous (2 + 1)- dimensional model one has f c (P) for 
(3 — > and 

/c( V-i + i>o g2 + (l). ,5,) 
Re-entrance does take place. 

Proof of Proposition \5.1l It is not difficult to prove that for the two-dimensional simple 
random walk 

G(y) = log(cosh(y) + 1) - log 2. (5.2) 

Indeed, first of all one has 

B[exp(yS^); S n / for n = 1, 2, . . . , N] < B[exp(yS^)] = Efexp^ 1 (5.3) 

= (1/2 + (1/2) cosh y) N 



for every N, so that 



G(y) < log(cosh(y) + 1) 



-log 2. 



(5.4) 
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To get the complementary lower bound, note that 

E[eK V (yS^);S n ^0 for n = 1, . . . , N] > E[exp{yS^); ^ for n = 1, . . . ,N}. (5.5) 

The quantity one is averaging in the right-hand side of (|5.5j) involves only {Sn^} n >o, 
which is just a (p, g)-walk with p = 1/2. Therefore, from the knowledge of G(y) for the 
(p, g)-walk (cf. Section ITT2*|) one obtains c(y) > log(l/2 + (1/2) cosh(y)) which concludes 
the proof of (jS2J). 

Next, it is known that K(2n) ~ c/(n(logn) 2 ) with c > while K(2n + 1) = 0. As 
a consequence (cf. (El Theorem 2.1]), f(/3) vanishes for (3 —* faster than any power of /3. 
This, together with the expansion G(y) = y 2 /4 + o(y 2 ) for y — > 0, implies that / c (/3) — ► 0. 
As for the large-/3 behavior, in analogy with the discussion in Section 12 a look at the 
dominant trajectories gives f(/3) = [3/2 + (logif(2))/2 + o(l//5), where now K{2) = 1/4. 
Since G(y) = y — 2 log 2 + o(l) for y — > oo, (|5.1j) follows from Proposition ll.il □ 

The disordered (2 + l)-dimensional model where uj n are IID symmetric random variables 
uj n = ±1 can also be analyzed in our framework, with the following result: 



Proposition 5.2. For the disordered (2 + 1)- dimensional pinning model, f c (P) 
(3 — > and 



Re-entrance is observed. 



21og2 + i^2-^logK(2j) 



for 



(5.6) 



The proof of Proposition (|5.2|) is essentially identical to that of Proposition 14.11 The 
reason why the factors 1/2 of (|4.1|) are replaced by 1/4 is just that our two-dimensional 
walk can touch the defect line only for n even, and therefore for large (3 dominant trajec- 
tories will touch half of the total positive charges, i.e., approximately A r /4 of them. The 
occurrence of re-entrance follows from a numerical evaluation of the constant multiplying 
1/(3 in ()5.6|) . which turns out to be equal to 1.2427. . .. 



Appendix A. Proofs and technical estimates 

A.l. Proof of Proposition ri.il and Lemma l4.3l Proof of Provosition M . 1\ Decompose 
the trajectory according to the location of the last visit of S to the origin before N 
(including N) getting thus by the Markov property of S 



N 



m=0 



exp (3^2u n l Sn =o ; S„ 



E[exp (f3fS N „ m ); S k ^O,k = l, 



m] . 
(A.l) 



The right-hand side has the form ^2 m=0 a m (uj)bN-m, with ao(iv) = bo = 1, a m (uj) = 
exp((F(/3) + o(l))m) for m — > oo, P(du;)-a.s. (the presence of the indicator function 
of the event S m = is irrelevant, see e.g. ,8s Remark 1.2]), and b m = exp((c(/3/) + 
o(l))m), again for m — > oo. At this point we remark that > max(aAr(o;), bj^) 

and therefore liminfjv_ >00 (l/./V) log Zjj* is bounded below by max (f((3), g(/?/)). For the 
opposite inequality we notice that for every e > there exists A(lv) such that a n (uj) < 
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A(u) exp ((f(/3) + e)n) for every n. Analogously, 6 n (o;) < B exp ((g(/3) + e)n) for some 
constant 1? and every n. Therefore 

AT 

< A(u)B exp (F((3)m + G((3f)(N - m) + eiV) , (A.2) 

m=0 

so that 

Z nI <{N + l)A(u)B exp (max (f(/3), g(/3/)) A^ + eN) , (A.3) 
for every A" and, since e > can be chosen arbitrarily small, the proof is complete. □ 



Proof of Lemma \4-3[ We set 



-. oo 

Q ■= -^2-n og ir(i), (A.4) 



and we separate the proof in lower and upper bound. 

For the lower bound we select the 5 trajectories hitting if and only if the charge is 
+1 on that site. This yields 

/AfivH \ / A^H \ 

Z Z * I II ex P (/3)if(^)j K I N- e A > ( A - 5 ) 

where l\ := inf{n > : u n = 1}, i^+i '■= inf{n > : w n+ ^fe e . = 1} and Nn{^) : = 
maxjfc : Ej=i ^' < N}. Note that 

is an IID sequence of geometric random variables 
of parameter 1/2. By taking the logarithm and dividing by A" both sides in (|A.5|) . in the 
limit as A" — > oo we get 

1 A/jv(u>) 
F(J3) > /^limsup-A^M+limsup- V logif(^), (A.6) 

AT-+oo A N-*oo A ~[ 

where we have used the fact that K (n - Zj) > K(N) and (log K(N))/N -» 

as A" — > 00. By the classical Renewal Theorem A/jv(u;)/A" actually converges P(dw)-a.s. 
to 1/E[£i] = 1/2 and, in turn, by the Strong Law of Large Numbers also the second 
superior limit in the right-hand side of (|A.6|) is an almost sure limit and it is equal to 
E[log A'(^i)]/E[^i], which coincides with Q. This concludes the proof of the lower bound. 

For the upper bound we define A^^ '■= {n '■ Wn = +1} H [1, N] and for every e > and 
any realization of lo we define the set of trajectories 

^N,u>,e ■= {S : |{1 < n < N : S n = 0} A A Nfi) \ < e\A N>UJ \} , (A.7) 

where A denotes the symmetric difference of sets. Of course, by the Law of Large Numbers 

\An,lj\/N ^^5° 1/2, P(du;)-a.s.. Therefore for the partition function restricted to fij\r W£ 
we have 

Z Z ( n N,^) ^ ex P (P\ A **>\0- - <0) N ~< H ^p (N^ (l - I)) , (A.8) 

for some Nq(u>) which is P(du)-a.s. finite. We can therefore focus on Z*f* (^N,u),e), which 
is bounded above by exp(/3|^4jv,tj|)P (^N,ui,e) and it is thus sufficient to show that 

lim sup lim sup — log P (fijv, w , e ) < Q, (A. 9) 

e\0 N-*oo A 
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to conclude. 

In order to establish (|A.9|) we introduce a coarse graining length L G N (L is sent to 
oo in the end, that is after N —* oo and e \ 0, so that, in particular, eL can be chosen 
arbitrarily small and we assume below that eL -C 1). We assume that N/L G N and 
we break {1, . . . , N} into ./V/L non-overlapping blocks {J5j}j=i jv/i, of length L. For 
every realization of the disorder u; we decompose the event £In,u,e into the disjoint union 
of the events ^In,u>,£,vi 

v G {0, 1} N/L , defined by the property that if S G fijv, u Eii then 
n {n : 5 n = 0} = Bj n Ajv ja; if and only if Uj = 0. In short, there is a mismatch in 
the block Bj with respect to the energetically optimal contact configuration if and only if 
Vj = 1 (the mismatch can be on a single site or on several sites). Note that there cannot be 
more than eN blocks containing a mismatch, that is if \v\ := \{i : V{ = 1}\ > eN the event 
^n,uj,e,v is empty. Since eL <C 1, only a small fraction of the blocks contains mismatches. 
We have therefore 

P(^W) = Y\ P^N^cv) < exp(c(eL)N/L) max P (O^,^) , (A.10) 

\v\<eN 

\v\<eN 

x\.0 

where c(x) — > comes from estimating the cardinality of {v : \v\ < eN}. 

We are left with estimating P (p,N,u,£,v) uniformly in v. For this we introduce the 
random variable Kj = inf{n G Bj : iv n = 1} (k^ = oo if w n = for every n £ Bj) and for 
every j such that Kj < oo the event 

E i,u : = i S : B 3 n{n: s n = o}n [ Kj ,N] = A NjU! n n [ Kj , N]} , (A.ll) 

that is simply the event that there is no mismatch in Bj from step Kj onward. If Hj = oo 
then is just the set of all possible polymer trajectories, without any restriction. 

Finally, we define 

, , fP (eP t ,\S k , = 6] if K* < oo, , 
Yj(io) := I V ' / (A.12) 

[ 1 otherwise. 

Note that {Yj}j is a sequence of IID random variables and that, by the Renewal Theorem 
(very much like in the proof of the lower bound: we are essentially evaluating the same 

quantity), (l/L) log Y\(uj) L ^°? P(dw)-a.s.. Therefore, for every 5 G (0,1) we can find 
Lq such that for L > Lq 

P(yi(w) < exp((Q + S)L)) > 1--. (A.13) 

Notice now that 

p {nN,u,,e,v) < p (n i:Wi =o4 lW ) . ( A - 14 ) 

and if we set k := \{i : vi = 0}| and {i : vi = 0} = . . . ({ii}i=i,...,k increasing) we 
have that if Ki k < oo 

oeQ = p (pU^L,. n [s. lk = o} n e%J < p J Y ik 



'\,r, n/-|.._ j = p (^n^^^ n |6' Kife = 0| n E^j < p ^n— r Jfc ^), 

(A.15) 

where we have used the Markov property of S and in the last step we have neglected 
the event {S Ki = 0}. With our definition of Yj(u), the factorization inequality in (|A. 15|) 
actually holds also for Ki k = oo, so that by iterating we obtain 

r \,:r, ) < II Y M- (A.16) 
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By putting equations (|A.13|) to HA,16|) together, by applying the Strong Law of Large 
Numbers and by exploiting the fact that more than a fraction (1 — eL) of the N/L blocks 
is free of mismatches, we obtain that 

P {n N ^ £ ,v) < exp ((1 - 5)(Q + S)L(1 - 2eL)^\ , (A.17) 

for N > Nq(lo), with Nq(lo) a random value that is P(dcj)-a.s. finite. 
Overall we have therefore established that P( dw)-a.s. 

limsupilogP^iv^e) < ^- + (l-S)(Q + 8)(l-2eL) £ -^(l-S)(Q + 5). (A.18) 

Since 5 can be chosen arbitrarily small the proof of the upper bound is complete. □ 

A. 2. Estimates in the periodic case. We prove here Theorem 13.11 The first result, 
that is the small temperature expansion in Q3.1|) . just comes from evaluating the ground 
state energy. We omit the details since they are substantially easier than the ones needed 
for the analogous result in the disordered set-up (cf. Lemma l4.3l and App. lA~Tj) . We point 
out also that such a result follows directly from the semi-explicit solution available for 
periodic models 011] that we now outline since we exploit it in order to establish the high 
temperature expansion (j3.2|) . 

The free energy of periodic models can be expressed by first introducing the (Abelian) 
group § := Z/ (t(lu)Z), that is {1, . . . , t(w)} with periodic boundaries. With abuse of 
notation an element a of § is going to be identified with a point in {1, . . . ,t(w)}, so by 
n £ a we mean n = kT(iv) + a for some k £ Z. For b > we set 

K a (b) := ^K(n)exp(-6n), (A.19) 

and in turn for a and 7 G § 

A^ a (b, f3) := K a _ y (b) exp (fiu a ) . (A.20) 

By the Perron-Frobenius theory on matrices with positive entries, the t(cj) x T(w)-matrix 
A(b, (3) has a maximal positive eigenvalue, often called spectral radius of A(b,/3), that we 
denote by \j(b,(3). By standard arguments one shows that A w (-,/3) is decreasing and 
smooth. In [HI E] it is shown that the (w-dependent) free energy is if X u (0, /3) < 1. If 
instead X w (0, (3) > 1 then there exists a unique solution b > to the equation X u (b, (3) = 1 
and such a 6 is precisely the free energy f((3). 

Let us now expand A(b, (3) for small values of b and (3: 

A %a (b, (3) = A^O, 0) + (3u; a + X -(3 2 + o{(3 2 )^j - ^| (b 1 ' 2 + o(6 1 / 2 )) , (A.21) 

where the last term follows from the Riemann sum approximation procedure 

(1 — exp(— bn)) 



£ (1 - exp(-6n)) K(n) = b^ 2 ( vW(^) + «(!)) b £ { - 



(A.22) 



and by the fact that the integral in the last term is equal to 2y/ir. We now use the 
fact that the maximal eigenvalue X(A + eB), A matrix with positive terms and e small, 
can be written up to 0(e 2 ) terms as A(^4) + eu ■ Bv, with u and v respectively right 
and left eigenvectors of A with eigenvalue X(A), normalized by setting ^2 a v a = 1 and 
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J2 a u a v a = 1. In our case A = A(0, 0) turns out to be bi-stochastic, so A(^4), before 
denoted A a) (0, 0), is equal to 1 and v a = 1/t(u;), as well as u a = 1 for every a. This leads 
to the expansion 



X(b, 0) = 1 + — - V uj n + -0 2 - ^b 1 ' 2 + r(b, 0), (A.23) 




where r(6,/3) = O{0 2 ) + o(V&) if En=l w ™ ^ and r ( 6 " 9 ) = o 0# 2 ) + otherwise. 
Therefore the existence of a (unique) solution 6 = f(0) to A w (6, /3) = 1 for small requires 

En=i > and 



p(/3) = J W^n=r ,r , — (A.24) 



8p 



+o(/3 4 ) if e;l>« = o 



Armed with these asymptotic behaviors, (j3.2|) follows from G x (y) = -\/2y(l + o(l))/p 

(y\o). 

What happens when Yln=i < is that A w (0,/?) is smaller than 1 for small (see 
(lA~23jn . And if we set /3 := sup{/3 : A w (0,/3) < 0} (note that < oo unless co n = — 1 for 
every n, as one can see from the large expansion) then one readily sees that f c {0) > 
for > 0q and f c (0) = otherwise. □ 
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